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Abstract 

The representation of solutions of Maxwell's equations as superpositions of scalar wave- 
lets with vector coefficients developed earlier is generalized to wavelets with polarization, 
which are matrix-valued. The construction proceeds in four stages: (1) A Hilbert space 
7i of solutions is considered, based on a conformally invariant inner product. (2) The 
analytic-signal transform is used to extend solutions from real space-time to a complex 
space-time domain T . The evaluation maps £ z , which send any solution F = B + zE to 
the values F(z) of its extension at points z G T, are bounded linear maps on TC. Their 
adjoints = £* are the electromagnetic wavelets. (3) The eight real parameters 
z = x + iy G T are given a complete physical interpretation: x = (x, i) G R 4 is 
interpreted as a space-time point about which \I/ z is focussed. The imaginary space-time 
vector y = (y, s) is time-like, i.e., |y| < \s\. The sign of s is interpreted as the helicity 
of the wavelet, while \s\ is its scale. The 3- vector v = y/s is the velocity of its center. 
Thus wavelets parameterized by the set of Euclidean points E = {(x, is)} (real space and 
imaginary time) have stationary centers, and wavelets with y ^ are Doppler-shifted 
versions of ones with stationary centers. All the wavelets can be obtained from a single 
"mother" by conformal transformations. (4) A resolution of unity is established in Tl, 
giving a representation of solutions as "atomic compositions" of wavelets parameterized 
by z G E. This yields a constructive method for generating solutions with initial data 
specified locally in space and by scale. Other representations, employing wavelets with 
moving centers, are obtained by applying conformal transformations to the stationary 
representation. This could be useful in the analysis of electromagnetic waves reflected 
or emitted by moving objects, such as radar signals. 



1. Introduction 

In this paper we further develop the wavelet formulation of classical electrodynamics 
begun in Kaiser [12-14]. There, it was shown that electromagnetic waves (solutions 
of Maxwell's equations) can be expressed as linear superpositions of spherical wavelets 
uniquely adapted to these equations. However, the wavelets in Refs. [12-14] were scalar- 
valued solutions of the wave equation rather than (vector- valued) solutions of Maxwell's 
equations. (Their coefficients were vector- valued, thus providing for polarization.) Al- 
though the scalar wavelets sufficed for the reconstruction of known electromagnetic 
waves, they could not be used for the construction of new ones according to given local 
data, precisely because of their scalar nature, which ignored the polarization degrees 
of freedom. (Their reproducing kernel was not the projection operator to the solu- 
tion space.) The wavelets constructed here remedy this deficiency, since they are true 
electromagnetic fields, parameterized by a complete set of physically relevant variables: 
Their point and time of localization, the velocity of their center, their scale and their 
helicity. 

This work is part of a general program whose main objective is to extend the el- 
ementary physical Gelds from real to complex space-time and interpret the imaginary 
parts of the space-time variables as control parameters for the wave number and fre- 
quency contents of the fields being analyzed. This approach has so far given a fully 
relativistic phase-space description of a variety of free field theories: Klein-Gordon, 
Dirac, and now Maxwell. Furthermore, the regularity resulting from the analyticity 
of the fields in the complex space-time domain promises to help resolve some of the 
fundamental difficulties plaguing these theories, related to their reliance on the concept 
of precise geometrical points with no substance or structure. For the extended fields, 
points in complex space-time have a natural interpretation as moving extended objects 
in real space-time which, in turn, act as elementary building blocks or "atoms" for the 
fields. In the case of field theories with positive mass (Klein-Gordon, Dirac), these 
atoms are relativistic coherent states, i.e., Gabor-like wavelets whose windows undergo 
scaling (Lorentz contractions) under Lorentz transformations. In the case of massless 
field theories such as electrodynamics, the atoms are space-time-scale wavelets trans- 
forming covariantly under the conformal group. In all cases, the atoms are uniquely 
determined by the field theory through covariance and analyticity. 

The main ideas of the above program were reported in [16], where some of the re- 
sults of the author's thesis (Phase-Space Approach to Relativistic Quantum Mechanics, 
Univ. of Toronto, 1977) were summarized. A key tool for extending general fields to 
complex space-time, the analytic-signal transform, was developed in [17] and further 
investigated in [11] . Some speculation on the application of these methods to elec- 
trodynamics was advanced in Kaiser and Streater [15] for the much simpler case of a 
two-dimensional space-time. 

In Section 2, Maxwell's equations are solved by Fourier methods from a viewpoint 
in which the concepts of helicity and polarization become very clear. A Hilbert space 
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Tt of solutions is constructed which was proved by L. Gross [8] to carry a unitary 
representation of the full invariance group of the equations, namely the conformal group 
of space-time. In Section 3 we review the analytic-signal transform, which extends any 
function f(x) from R n to C n . In general, the extended function f(z) is not analytic 
(there may not exist any analytic extension). But when the Fourier transform of / 
is supported on a double convex cone (as it is, for example, when / represents a free 
relativistic field such as an electromagnetic wave), then f(z) is analytic in a certain 
tube domain T in C n . In Section 4 we show that the analytic-signal transform, when 
applied to electromagnetic fields, uniquely determines a set of electromagnetic wavelets. 
A resolution of unity is derived which allows solutions to be expressed as superpositions 
of wavelets. In Section 5 we compute the reproducing kernel defined by the wavelets, 
which in turn gives the wavelets explicitly. In Section 6 we show how arbitrary solutions 
in 7i can be constructed from wavelets, with initial data specified locally and by scale. 
In Section 7, the wavelets are given a complete physical and geometric interpretation. 
In Section 8 we describe some generalizations, and a possible application. 

2. The Fourier Representation of Free Electromagnetic Fields 

An electromagnetic wave in free space (without sources or boundaries) is described by 
a pair of vector fields depending on the space-time variables x = (x, xq) (where x is the 
position and xq is the time), namely the electric field E(x) and the magnetic field B(x). 
These are subject to Maxwell's equations, 

V x E + d B = 0, V • E = 0, 

V x B — <9 E = 0, V • B = 0, ^ 

where V is the gradient with respect to the space variables and <9 is the time derivative. 
We have set the speed of light c = 1 for convenience. For the present dicussion, E and 
B may be assumed to be tempered distributions, so that (2.1) holds weakly. Later, the 
fields will be required to belong to a certain Hilbert space. Note that the equations are 
symmetric under the linear mapping defined by J : E i— > B, B h — E, and that J 2 is 
minus the identity map. Such a mapping on a vector space is called a complex structure, 
by analogy with multiplication by i in the complex plane. The combinations B ± i E 
diagonalize J, since J(B ± iE) = ±i(B±iE). They each map Maxwell's equations 
to a form in which the concepts of helicity and polarization become very simple. It 
will suffice to consider only F = B + i E, since the other combination is equivalent. 
Eqs. (2.1) then become 

d F = iV x F, V • F = 0. (2) 

Note that the first of these equations is an evolution equation (initial- value problem), 
while the second is a constraint on the initial values. Taking the divergence of the first 
equation shows that the constraint is conserved by time evolution. Note also that it is 
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the factor % in (2.2) (i.e., the complex structure!) which couples the dynamics of the 
electric field E to those of the magnetic field B. Eq. (2.2) implies 

-dlF = V x (V x F) = V(V • F) - V 2 F = -V 2 F (3) 

in Cartesian coordinates, hence the components of F become decoupled and each satis- 
fies the wave equation 

□ F = {-dl + V 2 )F = 0. (4) 
Since F is a tempered distribution, it has a Fourier expansion 

F(x) = (2tt)- 4 [ d 4 pe ip - x F(p), (5) 

where p = (p,po) G R 4 with p G R 3 as the spatial wave vector and po as the frequency. 
We use the Lorentz-invariant inner product p • x = PoXq — p • x. The wave equation 
(2.4) implies that p 2 F(p) = 0, where p 2 = p ■ p = pi — |p| 2 . If the distribution F has 
no essential support at the origin p = (i.e., no term proportional to S(p)), it must be 
supported on the nipped light cone 

C = {(p,po) G R 4 \{0} : p 2 = |p| 2 } = C+ U C-, (6) 

where ±po = |p| > in C±. Hence F has the form F(p) = 2nd(p 2 ) f (p), where f(p) 
is a (vector-valued) distribution on C or, equivalently, the pair of distributions on R 3 
given by f±(p) = f(p, ±|p|). For the moment, we assume that f± are (vector- valued) 
Schwartz test functions. Later the class of f 's will be enlarged by introducing an inner 
product and completing it to a Hilbert space, subject to a restriction related to our 
having "nipped" the light cone, which amounts, roughly, to f(0) = 0. Letting u> = |p|, 
we have 

8{p ) = 8{{po -oj)(po + oj)) = ^— i . (7) 

Hence (2.5) becomes 

F(x) = (27T)" 3 / |^e-^ x [e^f+Cp) + e- fat f_(p)] 

7r3 2 bo I - (8) 

dp e ip - x f(p), 

where dp = (27r) _3 (i 3 p/2 \p \ is a Lorentz-invariant measure on C. In order for (2.8) to 
give a solution of (2.2), f(p) must further satisfy the algebraic conditions 

ipof(p) = pxf(p), p-f(p) = (9) 
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for all peC, and the first of these equations suffices since it implies the second. Let 
v (p) = p/r>5 so that p G C if and only if |v| = 1. Define the operator T = T(p) on 
arbitrary functions g : C — > C 3 by 

(r g )(p) = -iv(p)xg(p), pec. (io) 

r(p) is represented by the Hermitian matrix 
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with matrix elements T mn (p) = ip 1 ^2 k=1 £ m nkPk , where e mn k is the totally antisym- 
metric tensor with £123 = 1. In terms of T(p), (2.9) becomes 

rf (p) = t{p). (12) 

Now for any g : C — > C 3 , 

T 2 g = -v x (v x g) = g - (v g )v, (13) 

so T(p) 2 is the orthogonal projection to the subspace of C 3 orthogonal to v(p), and it 
follows that 

r 3 g = r g . (14) 

The eigenvalues of r(p), for each p £ C, are therefore 1,0 and —1, and (2.12) states 
that f(p) is an eigenvector with eigenvalue 1. Since T(p) = —T(p), (2.12) implies that 
Tf = — f. A similar operator was defined and studied in much more detail by Moses 
[19] in connection with fluid mechanics as well as electrodynamics. 

Consider a single component of (2.8), i.e., the plane-wave solution 

F p (a;) = e ip - x f(p) = B p (x)+iE p (x), (15) 

with arbitrary but fixed p G C and f(p) 7^ 0. The electric and magnetic fields are 
obtained by taking the real and imaginary parts. Now T(p) f(p) = f(p) and T(p) f(p) = 
—f(p) imply 

T(p)F p (x) = F p (x), r(p)F7M=-FVM- (16) 

Since T(p)* = T(p), these eigenvectors of T(p) with eigenvalues 1 and —1 must be 
orthogonal: F p (x)*F p (x) = F p (x) - F p (x) = 0, where the asterisk denotes the Hermitian 
transpose. Taking real and imaginary parts, we get 

\B p (x)\ 2 = \F p (x)\ 2 , B p (x) ■ E p (x) = 0. (17) 
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The first equation shows that neither B p (x) nor E p (x) can vanish at any x (since 
f(p) 7^ 0). Furthermore, (2.16) implies that p x E p (x) = poTi p (x). Thus, for any x, 
{p, E p (x), B p (x)} is a right-handed orthogonal basis if po > (i.e., p G C+) and a 
left-handed orthogonal basis if po < (p E C_). Taking the real and imaginary parts 
of (2.15) and using f(p) = B p (0) + zE p (0), we have 

B p (x) = cos(p ■ x)Bp(O) - sin(p • x)E p (0), 
E p (x) = cos(p • x)Ep(O) + sin(p • x)B p (0). 

An observer at any fixed location x G R 3 sees these fields rotating, as a function of 
time, in the plane orthogonal to p. If p G C+, the rotation is that of a right-handed 
corkscrew, or helix, moving in the direction of p, whereas if p G C_, it is that of a left- 
handed corkscrew. Hence F p (x) is said to have positive helicity if p G C+ and negative 
helicity if p G C_. 

A general solution of the form (2.8) has positive helicity if f(p) is supported in 
C_|_ and negative helicity if f(p) is supported in C_. Other states of polarization, 
such as linear or elliptic, are obtained by mixing positive and negative helicities. The 
significance of the complex combination F(x) = B(x) +i E(x) therefore seems to be that 
in Fourier space, the sign of the frequency po gives the helicity of the solution! (Usually 
in signal analysis, the sign of the frequency is not given any physical interpretation, and 
negative frequency is regarded as a convenient mathematical artifact.) In other words, 
the combination B+z E "polarizes" the helicity, with positive and negative helicity states 
being represented in C+ and C_ , respectively. Had we used the opposite combination 
B — % E, C + and C_ would have parameterized the plane-wave solutions with opposite 
helicities. Nothing new seems to be gained by considering this alternative. (In fact, 
Maxwell's equations are invariant under the continuous group of duality rotations, of 
which the complex structure J mapping E to B and B to — E is a special case. In the 
complexified solution space, the combinations B ± % E form invariant subspaces with 
respect to the duality rotations. That gives the choice of B + zE an interpretation in 
terms of group representation theory.) 

In order to eliminate the constraint, we now proceed as follows: Let 



U(p)^ 1 -[F(p)+T^p)]. (19) 
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The established properties r* = T = T 3 imply that II* = II = II 2 and ITI = II, which 
proves that H(p) is the orthogonal projection to eigenvectors of F(p) with eigenvalue 
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1. Thus, to satisfy the constraint (2.12), we need only replace the constrained function 
f(p) in (2.8) by II(p)f(p), where now f(p) is unconstrained: 

F(x) = [ dp e ip ' x n(p) f (p) . (21) 
Jc 

Consequently, the mapping f i— > F is not one-to-one since II is a projection operator. In 
fact, f is closely related to the potentials for F, which consist of a real 3- vector potential 
A(x) and a real scalar potential Aq{x) such that 

B = V x A, E = -d A - VA . (22) 

The combination (A(x),Aq(x)) is called a "4-vector potential" for the field. We can 
assume without loss of generality that the potential satisfies the Lorentz condition 
V • A + doA = (Jackson [10]). Since A and A also satisfy the wave equation 
(2.4), they have Fourier representations similar to (2.8): 

A(x)= [ dpe ip - x a{p), A (x)= [ dp e ip ' x a (p). (23) 
Jc Jc 

The Lorentz condition means that poao(p) = p ■ a(p), or ao(p) = v • a(p), so ao is 
determined by a. Eqs. (2.22) will be satisfied provided that the Fourier representatives 
e(p), h(p) of E, B satisfy 

b = -ipxa = p ra, e = — zp a + i pao = — ipo T 2 a. (24) 
Hence F = B + % E is represented in Fourier space by 

b(p) + z e(p) = po [r(p) + T(p) 2 ] a(p) = 2p n(p) a(p). (25) 

This shows that we can interpret the unconstrained function f(p) in (2.21) as being 
directly related to the 3-vector potential by 

f(p) = 2p a(p), (26) 

modulo terms annihilated by II (p), which correspond to eigenvalues —1 and of r(p). 
Seen in this light, the non-uniqueness of f in (2.21) is an expression of gauge freedom 
in the B + i E representation, as seen from Fourier space. In the space-time domain, 
(B, E) are the components of a 2-form F in R 4 and (A, Aq) are the components of a 
1-form A. Then Eqs. (2.1) become dF = and 5F = (where 5 is the divergence with 
respect to the Lorentzian inner product), Eqs. (2.22) become unified as F = dA, the 
Lorentz condition reads 5 A = and the gauge freedom corresponds to the invariance of 
F under A — > A + dx, where xi x ) is a scalar solution of the wave equation. 

Maxwell's equations are invariant under a large group of space-time transforma- 
tions. Such transformations produce new solutions from known ones by acting on the 
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underlying space-time variables (possibly with a multiplier to rotate or scale the vector 
fields). Some trivial examples are space and time translations: Obviously, a translated 
version of a solution is again a solution, since the equations have constant coefficients. 
Similarly, a rotated version of a solution is a solution. A less obvious example is Lorentz 
transformations, which are interpreted as transforming to a uniformly moving reference 
frame in space-time. (In fact, it was in the study of the Lorentz invariance of Maxwell's 
equations that the Special Theory of Relativity originated; see Einstein et al. [6].) The 
scale transformations x — > ax, a ^ 0, also map solutions to solutions, since Maxwell's 
equations are homogeneous in the space-time variables. Finally, the equations are in- 
variant under "special conformal transformations" (Bateman [2], Cunnigham [4]), which 
can be interpreted as transforming to a uniformly accelerating reference frame (Page 
[20]; Hill [9]). Altogether, these transformations form a 15-dimenional Lie group called 
the conformal group, which is locally isomorphic to 577(2, 2) and is here denoted by C. 
Whereas wavelets in one dimension are related to one another by translations and scal- 
ings, electromagnetic wavelets will be seen to be related by conformal transformations, 
which include translations and scalings. (A study of the action of 577 (2, 2) on solutions 
of Maxwell's equations has been made by Riihl [21].) 

To construct the machinery of wavelet analysis, we introduce a Hilbert space struc- 
ture on the solutions. It is important to choose the inner product to be invariant under 
the largest possible group of symmetries, since this allows the largest set of solutions in 
H to be generated by unitary transformations from any one known solution. (In quan- 
tum mechanics, invariance of the inner product is also an expression of the fundamental 
invariance of the laws of nature with respect to the symmetries in question.) Let f(p) 
satisfy (2.12), and let a(p) be a vector potential for f satisfying the Lorentz condition, 
so that the scalar potential is determined by a (p) = v ■ a(p). By (2.25), 

|f(p)| 2 = Ap 2 |n(p) a(p)| 2 = Ap 2 aT^) • U(p) a(p) 

= 2pl aT^) • T{p) a(p) + 2p 2 aT^) • T(p) 2 a(p). 

The first term is 



-2ipl a(p) • (v x a(p)) = 2ip 2 v • (a(p) x a(p)), 



(28) 



which cancels its counterpart with p 
8l(p). Thus 



-p on account of the reality condition a(— p) = 



dp 
c Po 



f (p) 1 2 = 2 / dp a(p) • [a(p) - v(v • a(p))] 



dp \&(p)\ 2 - \ao(p)\ 



= 2 / 

Jc 



(29) 



The integrand in the last expression is the negative of the Lorentz-square of the 4- 
potential (a(p), ao(p)). Consequently, the integral can be shown to be invariant under 
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Lorentz transformations. (Note that |a| 2 — |a | 2 > 0, vanishing only when a(p) is a 
multiple of p, in which case f = 0. This corresponds to "longitudinal polarization.") 
Hence (2.29) defines a norm on solutions which is invariant under Lorentz transforma- 
tions as well as space-time translations. In fact, the norm (2.29) is uniquely determined, 
up to a constant factor, by the requirement that it be so invariant. Moreover, Gross [8] 
has shown it to be invariant under the full confer mal group C. Again we eliminate the 
constraint by replacing f(p) with H(p) f(f>). Thus, let 7i be the set of all solutions F(x) 
defined by (2.21) with f : C — > C 3 square-integrable in the sense that 

||F|| 2 = f ^|II(p)f(p)| 2 = (27r)- 3 I ^|II(p)f(p)| 2 <oc, (30) 
Jc Po Jc 2 IpI 

H is a Hilbert space under the inner product obtained by applying the polarization 
identity to (2.30) and using (IIf)*IIg = PITTIg = f*IIg: 

(F,G)= / %f( P ru(p) g (p). (31) 
Jc Po 

H will be our main arena for developing the wavelet analysis and synthesis of solutions. 
Note that when (2.12) holds and f±(p) = f(p, ±|p|) are Schwartz test function as we 
assumed earlier, then 

f(0) = f±(O) = O (32) 

must hold in order that (2.30) be satisfied. However, now that we have our Hilbert 
structure, we complete to the larger class of all (generalized) functions f satisfying 
(2.30). 

To show the invariance of (2.30) under conformal transformations, Gross derived 
an equivalent norm expressed directly in terms of the values of the fields in space at any 
particular time xq = t: 

l|F||Lo SS = h [ F(x, t)*F(y, t). (33) 

71 Jrs i x - y| 

The right-hand side is independent of t due to the invariance of Maxwell's equations 
under time translations (which is, in turn, related to the conservation of energy). A 
disadvantage of the expression (2.33) is that it is non-local, since it uses the values of 
the field simultaneously at the space points x and y. In fact, it is known that no local 
expression for the inner product can exist in terms of the Geld values in (real) space- 
time R 4 (Bargmann and Wigner [1]). In Section 4, we derive an alternate expression for 
the inner product directly in terms of the values of the electromagnetic fields, extended 
analytically to complex space-time. This expression is "local" in the space-scale domain 
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(rather than in space alone). But first we must introduce the tool which implements 
the extension to complex space-time. 

3. The Analytic-Signal Transform 

Given a vector function F : R n — > C m , we define its analytic-signal transform as the 
function F : C n — > C m given by the following line integral in R n : 

1 f°° (It 

F(x + iy) = — :F(x + ry). (3.1) 

7U J_ oc t-i 

This transform was introduced and studied in Kaiser [11, 17], where it was shown to be 
related to the Fourier transform F(p) of F(x) by 

F(x + iy) = (2yr)" n I d n p29{p ■ y) e ip < x+iy) F(p) . (3.2) 
Jr.™ 

Here 9 is the unit step function, defined by 9{u) = if u < 0, 9(0) = |, 9{u) = 1 
if u > 0. For (3.2) to make sense, it suffices that F(p) be absolutely integrable, since 
\9(p-y) e ip '( x+w >\ < 1. For concreteness, we assume for the time being that F E L 1 (R 4 ) 
and use (3.2) to define F, viewing (3.1) as "motivation." (A study of (3.1) in the context 
of distribution theory is currently being undertaken by T. Takiguchi [24].) Note that 
setting y = on the right gives the inverse Fourier transform of F(p), so that formally 
we have F(x) = F(x) and F is an extension of F from R n to C n . (This is made more 
precise below, in (3.7).) Of course, this extension is usually not analytic, since in general 
there exists no analytic extension. However, when R n is space-time and F represents 
a free physical field such as an electromagnetic field (m = 3), a Klein-Gordon field 
(m = 1) or a Dirac field (m = 4), then F(p) = F(p,p ) vanishes outside the solid light 
cone 

V = {(p,po) e R 4 : P 2 = Vl ~ |p| 2 > 0, po ^ 0} = V+ U V-, (3.3) 

where ±p > |p| > in V±. (In the electromagnetic case, for example, F is supported on 
the boundary C of V as a consequence of the wave equation (2.4).) Hence the integral 
in (3.2) extends only over V. Formally, the obstacle to the analyticity of F(z) in (3.2) 
is the factor 9(p- y) (which is necessary, generally, to ensure that the integral converges 
in the region of Fourier space where e~ p ' y > 1). However, when F is supported in V, 
that obstacle can be removed as follows: Suppose y is such that p ■ y > for all p G V + 
and p ■ y > for all p E V_. (This means that the hyperplane Lorentz-orthogonal to y 
separates V + and V_!) Then 9(p ■ y) = 1 for all p E V + and 9(p ■ y) = for all p E V_, 
hence the integral now extends only over V+, and the obstracting factor is identically =1 
in that cone. Furthermore, the extra factor e~ p ' y coming from the analytic continuation 
of the Fourier kernel provides exponential damping, which leads to the analyticity of F 
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at (x + iy), for all x. Similarly, if y is such that p ■ y > for all p E V_ and p ■ y < 
for all p e V+, then F is again analytic at x + for all x. The above sets of imaginary 
space-time points y are, by definition, the dual cones V± of V± (Stein and Weiss [23]), 

V± = {y = (y, y ) eK 4 : p-y = p y - p • y > for all p e V±}. (3.4) 

and V'_ are the open future light cone and the open past light cone in space-time 
(as opposed to Fourier space, where V± live). The union V = U V'_ will be called 
the dual cone of V = V + U V_ . Explicitly, 

v£ = {(y, y ) : ±yo > |y|}, ^ = {y e R 4 : y 2 = yl - |y| 2 > 0}. (3.5) 

The argument used above to motivate the definition of V± can be made precise, leading 
to the conclusion that F(z) is indeed analytic in 

T = {z = x + iyeC 4 :yeV'} = T + UT_, (3.6) 

where T± is the set of z : s with y e V±. The fact that F and V are almost identical is 
due to our arbitrary choice c = 1 for the speed of light. V is actually "reciprocal" to 
V: As c increases, V narrows and V widens. 

A general function F(x) with suppF(p) C V therefore becomes "polarized" when 
extended to T: The positive-frequency part (p G V+) determines F(z) in T + , and 
the negative-frequency part determines it in 71. The positive and negative frequency 
components mix on the common boundary R 4 of 7+ and 71. (The boundaries of 7± 
are seven-dimensional, but their intersection is R 4 .) If only real vector functions F 
are considered, then the positive and negative frequency parts are "coupled" by the 

reality condition F(— p) = F(p), with the corresponding relation F(z) = F(z) in T. If 
F is allowed to be complex- valued, its positive- and negative- frequency parts become 
independent. In the case of electrodynamics, we saw in the last section that they 
correspond to the positive-helicity and negative-helicity parts of an electromagnetic 
wave in the B + zE representation. The separation of helicities into C + and C_ in 
Fourier space is translated, by the analytic-signal transform, to their separation into T + 
and 71. 

From a mathematical point of view, it suffices for F(p) to be supported in any 
double cone of the forrn V+ U V- , where V± are convex cones intersecting only at their 
common vertex. Then F is analytic at x + iy whenever the hyperplane orthogonal to y 
separates V + and V_, which again means that y belongs to the dual V of V, defined as in 
(3.5). The name "analytic-signal transform" derives from the fact that when n = m = 1 
and / is real- valued, then f(x + iy) coincides with the "analytic signal" of / for y > 0, 
as first defined by D. Gabor [7] in his famous paper on communication theory. In fact, if 
n = 1, then f(x + iy) is simultaneously the analytic extension of the positive- frequency 
part of / to the complex upper half-plane and of the negative-frequency part of / to 
the lower half-plane. (These two half-planes now play the roles of T + and 71.) 
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As already mentioned, F(x) = F(x) formally, i.e., by setting y = in (3.2). More 
precisely, F is the boundary value of F in the sense that for any y E R n \{0}, 



lim F(x + iey) + F(x — iey) = 2F(x) a.e. 



(3.7) 



On the other hand, the "jump" of F across R n is 



lim F(x + iey) 



F(x-iey) = 2iH y F(x) a.e., 



(3.8) 



where 




F(x - ry) 



(3.9) 



is the multidimensional Hilbert transform of F in the direction ofy^O (Stein [22]) and 
PV denotes the Cauchy principal value. 

4. The Electromagnetic Wavelets \I/ z 

We are now ready to pursue our main theme, the construction of the electromagnetic 
wavelets and their resolution of unity. (For general background on wavelet theory, the 
reader may consult Chui [3], Daubechies [5], Kaiser [18] and the references therein.) 
Consider the extension of the electromagnetic field F(x) to the tube domain T defined 
in (3.5) and (3.6). Combining (2.21) and (3.2), we obtain 



As earlier, assume that f±(p) = f(p, ±|p|) are vector-valued Schwartz test functions, 
to begin with. Then F is analytic in T. Fix an arbitrary z = x + iy E T (i.e., 
y 2 = y 2 — |y| 2 > 0) and consider the linear operator S z : H — > C 3 defined by £ Z F = F(z). 
This is an evaluation map which, when applied to the field F, gives the value of its 
extension at the complex space-time point z. Because of the analyticity of F, S z turns 
out to be bounded, as will be seen later. (It becomes unbounded as y 2 — > 0.) We now 
define the electromagnetic wavelets as the adjoint operators ^ z = S* : C 3 — > H. To find 
these explicitly, choose any orthonormal basis 111,112,113 of C 3 and let *& Z) k = ^ z^k £ 
7i, k = 1,2,3. This gives three solutions of Maxwell's equations, all of which will be 
wavelets "at" z. ^ z is a matrix-valued solution of maxwell's equations, obtained by 
putting the three (column) vector solutions *& Zy k together. It will be convenient to use 
the following "star notation," intoduced in Kaiser [17]: For any F E H, let F* : H — > C 
denote the linear functional obtained by taking inner products with F: 



F(x + iy)= [ dp 2e(p-y)e ip - ix+iy) U(p){(p) 



Jc 



(4.1) 



F*G = (F,G), 



g e n. 



(4.2) 
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F* is not to be confused with the Hermitian transpose F(x)* of F(x) G C 3 . Then the 
k-th component of F(z) with respect to the basis {u fc } is 

F k (z) = utF(z) = u* k £ z F = utVlF = (* ,u fc )*F = ( * 2fc ,F). (4.3) 

By (4.1), 

u* k F(z) = [ %2p 2 9(p ■ y) a** u* U(p) f (p) , (4.4) 
which shows that ^ Zjk is given in the Fourier domain by 

^ >fc (p) = 2p 2 <9(p . y ) e -*P-* n(p) u fc . (4.5) 

Note that each i/) zk (p) satisfies the constraint since T(p)U(p) = II (p). The matrix- 
valued wavelet *Sf z in the Fourier domain is therefore 

^(p) = 2^0(p.y)e-*-*II(p). (4.6) 
In the space-time domain we have (using H(p) i/> z (p) = ij) z {p)) 

* z (x') = [ dp e^'^M = [ dp 2p 2 9(p-y)e^'-^U(p). (4.7) 
Jc Jc 

Now that we have the wavelets, we want to make them into a "basis" that can 
be used to decompose and compose arbitrary solutions. This will be accomplished by 
constructing a "resolution of unity" in terms of the wavelets. To this end, we derive an 
expression for the inner product in H directly in terms of the values F(z) of the extended 
fields. To begin with, it will suffice to consider the values of F only at Euclidean space- 
time points, i.e., at points with a imaginary time coordinate zq = is and real space 
coordinates z = x. In order for z to belong to T, it is only necessary to have s ^ 0. 
We denote the set of all such points by E. The name "Euclidean" stems from the fact 
that at such points, the negative of the indefinite Lorentzian metric restricts to the 
positive-definite Euclidean metric on E: —z 2 = —(is) 2 + |x| 2 = s 2 + |x| 2 . Later, x 
will be interpreted as the center of the wavelets ^ Zjk , and s as their helicity and scale. 
Using (4.1) and letting ui = |p| = |po|> we have 

F(x,is) = / dp26(p s)e- poS - ip - x Il(p)f(p) 
Jc 

= 2 [ dp e" ip ' x [6(us) e- ws II(p, u) f (p, to) + 6(-us) e us II(p, -ui) f (p, -u)] 
u}- x 0(us) e- ws n(p, u) f (p, u) + uj-^i-us) e ws TL(p, -u) f (p, (x), 

(4.8) 
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where v denotes the inverse Fourier transform with respect to p. Hence by Plancherel's 
formula, 



/ d 3 x|F(x,z 

</R3 

d 3 p 



'R3 

r 

-2u>s Itt/„ , .\ p/„ , . \ 1 2 i /)/ , ,„\ J2los |tt/„ , .\ p/„ , .\|2 



r d 6 

7r3 (2tt) 



/ R 3 (2tt)V 

where we used 6(u) 2 = 6(u) and 9(u) 9{—u) = for u ^ 0. Thus 



0(tus) e~ Zujs |n(p, u) f (p, to)\ z + 6{-us) e Zujs |n(p, -w) f (p, -u) 

(4.9) 



/ rf 3 x^|F(x,z S )| 2 = f ^ \ \U(p, u) f (p, c^)| 2 + |n(p, -a;) f (p, -«,) | 2 
Je Jrs 2(27t) a; J L 

= / ^-|i%)f( P )| 2 = / ^f(p)*n(p)f( P ) = ||F|| 2 , 

./c Po Jc Po 

(4.10) 

since n*n = FE 2 = FL Let H be the set of all analytic-signal transforms F of solutions 
FeH. For F, G G U, write 

(F, G) = I d 3 x ds F(z)* G(z). (4.11) 
Then (4.10) leads immediately to the following result. 

Theorem 1. H is a Hilbert space under the inner product (4.11), and the map F i— > F 
is unitary from 7i onto H. 

Proof: By the polarization identity, (4.10) implies 

(F,G) = (F,G), (4.12) 
so the map is an isometry. It is obviously surjective, by the definition of 7i. I 

With the "star notation" introduced earlier, the Hermitian transpose F(z)* : C 3 — > C 
of the 'column vector' F(z) G C 3 is the composition 

F(z)* = (**F)* = F** 2 , (4.13) 
where F* : Ji — > C denotes the linear functional (4.2). Hence the integrand in (4.11) is 

F(z)* G(z) = (** F)* ** G = F* * 2 G, (4.14) 
where \I/ z : H — > 7i is the composition of and \I/ 2 , and (4.11) reads 
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/ rf 3 x ds F* \I> Z G = F* G, F, G G W. 



(4.15) 



Theorem 2. 

(a) The wavelets \I/ 2 with z G E give the following resolution of the identity I in H: 

[ d 3 x ds * 2 *: =1, (4.16) 

where the equality holds in the weak topology ofH, i.e., (4.15) is satisfied. 

(b) Every solution F G H can be written as a superposition of the wavelets * z with 
z = (x, is) G E, according to 

F = f d 3 x ds xis ** is F = / rf 3 x ds xis F(x, is), (4.17) 
Je Je 

i.e., 

F(x') = / d 3 x * Xils (s)F(x,«s) a.e. (4.18) 
Je 

(4.17) holds weakly in TL (i.e., the inner products of both sides with any member ofTC 
are equal). However, for the extended fields, we have 

F(z') = V* Z ,F = [ d 3 x ds * XJS F(x, is) (4.19) 
Je 

pointwise for all z' G T. 

Proof: Only the pointwise convergence in (4.19) remains to be shown. This follows from 
the boundedness of which will be proved in Section 5. I 

The pointwise equality fails, in general, for the boundary values F(x) because the 
evaluation maps (or, equivalently, their adjoints * 2 ) become unbounded as y — > 0. This 
will be seen in the next section. 

The opposite composition \f r *,\f r z : C 3 — > C 3 is a matrix- valued function on T x T: 



K(z' | z) = **, * z = / ^4p<*% V)%-y)e^ (2 '- z_) I%) 2 

= 4 /dp rftfCp-^^-yJe^-^nCp). 

./c 



(4.20) 
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Eq. (4.19) shows that K(Y \z) is a reproducing kernel for the Hilbert space H; see 
Kaiser [11] for background and references. The boundary value of K(V \z) as y' — > 
is, according to (3.7) and (4.7), given by 

K(x' \z) = - lim [K(x' + isy' \ z) + K(x' - iey' | z)} = V z (x'). (4.21) 

Hence, to find the wavelets explicitly, we must compute their reproducing kernel. This 
is done in the next section. 

The meaning of the index k in *& Z) k deserves to be examined. Since Tl(p) is the 
orthogonal projection to the eigenspace of T(p) with the nondegenerate eigenvalue 1, all 
the columns (as well as the rows) of Tl(p) are all multiples of one another. But the coef- 
ficients are p-dependent, and the algebraic linear dependence in Fourier space translates 
to a differential equation in space-time, relating the different wavelets *f? z ,k- For the 
columns, this differential equation is just Maxwell's vector equation (2.2). (Recall that 
the scalar equation is then implied by the wave equation.) Since II (p) is Hermitian, the 
same argument goes for the rows. Explicitly, 

r(p)*,(p) = + z ip) = 4> z (p)t( p ). (4.22) 

When multiplied through by po and transformed to space-time, these read 

V x V z (x') = -id' V z (x') = Vr z (x') x V 7 , (4.23) 

where d' Q denotes the partial with respect to x' Q , V' the gradient with respect to x', and 

V' indicates that V' acts to the left, i.e., on the column index. This states that not only 
the columns, but also the rows of *& z are solutions of Maxwell's equations. The three 
wavelets ^ Zy k are thus coupled. Note also that since * z {x') = ^ z - x '(0), Eq. (4.23) can 
be rewritten as 

V x = -id * z = * 2 x V, (4.24) 

where <9 and V are the corresponding operators with respect to the labels x = Re z 
and x = Re z . 

We will see in Section 7 that the reconstruction of F(x') from F(x, is) can be 
obtained by a much simpler method than (4.18), using only a single scalar wavelet 
^x.,is(x') instead of the matrix wavelet ^ x ,is(x') (or three vector wavelets ^ x ,is,k(x')). 
However, that presumes that we already know F(x, is), and without this knowledge 
the reconstruction becomes meaningless, since no new solutions can be obtained this 
way. The use of matrix wavelets will be necessary in order to give a generalization of 
(4.17)-(4.19), where F(x, is) can be replaced with an unconstrained coefficient function. 
In other words, we need matrix wavelets in space-time for exactly the same reason that 
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II(p) was needed in Fourier space (Eq. (2.21)): To eliminate the constraints in the 
coefficient function. 



5. The Reproducing Kernel 

In order to obtain detailed information on the wavelets, we compute the reproducing 
kernel (4.20) explicitly. Note, first of all, that if y' • y < (i.e., z' G 7+ and z G 71 or 
z' £ T_ and 2 G 7+), then K(z' | z) = since p • j/' and p • y have opposite signs for 
all p G C. Hence it suffices to compute the kernel for z' and z in the same half of T. 
Furthermore, H(—p) = H(p) since T(—p) = T(p). Hence, letting z' — > z' and 2; — > z in 
(4.20) gives 



K(z / |z)=4 / dp p 2 e(-p-y')e(-p-y)e ip - (E '- z) U(p) 
Jc 

= 4 / dp p 2 o e{p-y')0{p-y)e ip - {z -^U(p) = K(z\z'), 
Jc 



(5.1) 



where the last equality is obtained by letting p — > —p. Thus it suffices to compute the 
kernel for z' , z G T + . In this case, 

K(z' I z) = 4 / dp e i P -(«'-2) n ( p ) = L ^ _ ^ ( 5 2 ) 

is analytic in «; = z' - z G T + . It can be shown that L(zy) with 1/ G K[ uniquely 
determines L(w) for all w G 7+ by analytic continuation, hence it suffices to compute 
only ~L(iy) for t/G V{. Now the matrix elements of 2pl U(p) are given by (2.20): 



2plU mn (p) = 5 mn pl -pmPn + i^2e mn kPoPk- (5.3) 

fc=l 

To compute L(zy), it is useful to write the coordinates of y in contravariant form: 

2/° = 2/o, 2/ m = -2/m ( m = 1,2,3), so that p ■ y = E)LoM M - Letting d M denote the 
partial derivative with respect to y^ for \i = 0, 1, 2, 3, we have 

2/ dp PuPve-™ = d tl d v S(y) = S IM1 ,(y), ^ = 0,1,2,3, (5.4) 



where 



5(2/) = 2 / dp e-P'f, yeVi. (5.5) 



Thus (5.3) and (5.4) give the matrix elements of L(iy) as 
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3 

L mn {iy) = S mn Soo(y) - S mn {y) + i ^e mn fc Spkjy), m,n = 1,2,3. (5.6) 

fc=i 

It only remains to compute S'(y). For this, we use the fact that S(y) is invariant under 
Lorentz transformations, since p ■ y and dp are invariant and C + is a homogeneous 
space for the proper Lorentz group. Since y G Vi, there exists a Lorentz transformation 
mapping y to (0, A), where X(y) = (y 2 , — ly] 2 ) 1 / 2 > 0. The invariance of S implies that 
S(y) = S(0, A). Letting u> = |p| again, we thus have 

Siy) = 2 L W" AlPl = ^1°°"^ = ^ (6 ' 7) 

Taking partials with respect to y M and y" gives 

S„{y) = Ay * Vv ~* lv>? , ^ = 0,1,2,3, (5.8) 
where g^ u = diag(l, —1, —1, —1) is the Lorentz metric. It follows that 



*2 r ^ 

L mn (iy) = S mn (y 2 +y\ + yl + yf ) - 2y m y n + 2i e mnk y y k 



(5.9) 



k=i 



To compute L(io) for w e T + , we need only replace y with —iw. This gives 



2 r d n 

£mn(«0 = o fT 5 mn (wo + ty? + w\ + w|) - 2w m w n + 2i V S mnk W W k , (5.10) 



fc=l 



where = («; • w) 3 . The full kernel is obtained by setting w = z' — z and multiplying 
by 9(y' ■ y), which ensures that it vanishes when z' and z are in opposite halves of T: 

K(z'\z)=e(y'-y)L(z'-z), z',zeT. (5.11) 

In Section 4 we stated that due to the analyticity of F(z), the evaluation maps S z 
(and with them, the wavelets \I/ z = £*) are bounded, and that they become unbounded 
as z = x + iy approaches the boundary of T, i.e., y 2 — > 0. This can now be verified by 
examining K(z \ z) = ** * 2 . By (5.11), 

K(^) = % 2 )L(2zy)=L(2zy) (5.12) 

for all z G T, since y 2 > in V'. Eq. (5.9) shows that \I/ 2 is indeed bounded when 
z G T diverges as y 2 — > 0. For example, if y = (0, s) (which can always be arranged by 
applying a Lorentz transformation), then 
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(5.13) 



where I is the identity matrix in C 3 . 

6. Atomic Composition of Electromagnetic Waves 

The reproducing kernel computed in the last section can be used to construct elec- 
tromagnetic waves according to local specifications, rather than merely to reconstruct 
known solutions from their analytic-signal transforms on E. This is especially interest- 
ing because the Fourier method for constructing solutions (Section 2) uses plane waves 
and is therefore completely unsuitable to deal with questions involving local properties 
of the fields. It will be shown in Section 7 that the wavelets i Sf x+ i y (x') are localized 
solutions of Maxwell's equations, at the "initial" time x' = x . Hence we call the 
composition of waves from wavelets "atomic." 

Suppose F is the analytic-signal transform of a solution F e H of Maxwell's equa- 
tions. Then according to (4.10), 



Let C 2 (E) be the set of all measurable functions : E — > C 3 for which the above 
integral converges. C 2 (E) is a Hilbert space under the obvious inner product, obtained 
from (6.1) by polarization. (In fact, we could identify E with R 4 and C 2 (E) with 
L 2 (R 4 ) since the set R 4 \-E = {(x,0):x£ R 3 } has zero measure in R 4 . But this could 
cause confusion between the Euclidean region E and real spacetime R 4 .) Define the 
map Re '■ 7~t — > C 2 (E) by 



That is, Re F is the restriction F \ e to E of the analytic-signal transform F of F. Then 
(6.1) implies that the range W of Re is a closed subspace of £ 2 (E), and Re maps Ti 
isometrically onto W. (In the Physics literature, an operator which transforms fields in 
real space-time to their counterparts in Euclidean space-time is called a Wick rotation.) 
The following theorem characterizes the range of Re and gives the adjoint R* E . 




(6.1) 



(i? E F)(x,z S ) = *; js F = F(x,z S ). 



(6.2) 



19 



Theorem 3. 

(a) The range of Re is the set W of all # G C 2 (E) satisfying the "consistency condition" 

$(V)= / rf 3 x dsK(z' |x,-is)$(x,is), (6.3) 

J E 

pointwise in z' G E. 

(b) The adjoint operator R* E : C 2 (E) -> 7i is given by 

i2^*= / d 3 xds* x , ls #(x,is), (6.4) 
where the integral converges weakly in Ti. 

Proof: If $ G W, then <fr(x, is) = F(x, is) for some F G 7i, and (6.3) reduces to (4.19), 
which holds pointwise in z' G E. On the other hand, given a function $ G C 2 (E) which 
satisfies (6.3), let F denote the right-hand side of (6.4). Then for any G G H, 

G*F= I d 3 x dsG(x,is)*$(x,is) = (R E G,&) C 2, (6.5) 
Je 

where we have used G** X)is = (** is G)* = G(x, is)*. Hence the integral in (6.4) 
converges weakly in Ti. The transform of F under Re is 

(R E F)(z') = ** z ,F= j rf 3 xrfs*:,* x , ls #(x,is) 

Je (6.6) 

= / d 3 x dsK(z' | x, -is) *(x, is) = &(z'), 
Je 

by (6.3). Hence <fr G W as claimed, proving (a). Eq. (6.5) states that {G,F)w = 
{ReG,&)c2. That shows that F = R E &, proving (b). I 

Eq. (6.4) constructs a solution R E <& G Ti from a coefficient function $ G C 2 (E). 
When $ is actually the transform -R^F of a solution F G Ti, then <fr(x, is) = * x is F 
and 

R* E * = [ rf 3 x ds * X)JS ** is F = F, (6.7) 



E 



by (4.17). Thus R* E R E = I, the identity in Ti. (This is equivalent to (6.1).) We now 
examine the opposite composition. 

Theorem 4. The orthogonal projection to W in C 2 (E) is the composition 
P = R E R* E : £ 2 ( E ) -»• ^ 2 (^)> which is given by 

(P^)(z')= [ d 3 xdsK(z'\x,-is)®(x,is). (6.8) 
Je 
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Proof: By (6.4), 



(R E R E &)(z') = ** z ,R E & = [ rf 3 x^*:,* x ,, s #(x,z S ) = (P#)(/), (6.9) 

J E 

since \f , */\E , X) i s = K(z', x, — is). Hence ReR e = P- This also shows that P* = P. 
Furthermore, R* e Re = I implies that P 2 = P, hence P is indeed the orthogonal 
projection to its range. It only remains to show that the range of P is W. If $ = 
Re F G W, then R E R* E <& = ReR* e Re F = i?£ F = Conversely, any function in the 
range of P has the form $ = R E R* E ® for some G C 2 (E), hence «fr = PeF where 
F = R* E e G 7Y. I 

When the coefficient function <fr in (6.4) is the transform of an actual solution, 
then R E reconstructs that solution. However, this process does not appear to be too 
interesting, since we must have a complete knowledge of F to compute F(x, is). For 
example, to compute F(x, is) by (3.1), we must know F(x, t) for all x and all t. Hence, no 
"initial-value problem" is solved by (6.4) when applied to # G W. However, the option 
of applying (6.4) to arbitrary $ G jC 2 (E) is a very attractive one, since it is guaranteed 
to produce a solution without any assumptions on $ other than square-integrability. 
It is appropriate to call R E the construction operator associated with the resolution of 
unity (4.16). It can be used to construct solutions in H from unconstrained functions 
# G jC 2 (E). In fact, it is interesting to compare the wavelet construction formula 

F(x')= / rf 3 x rfs* X)JS (x')*(x,zs) (6.10) 

J E 

directly with its Fourier counterpart (2.21): 

F(x')= [ dpe ip - x 'u(p)i(p). (6.11) 
Jc 

In both cases, the coefficient functions (<& and f) are unconstrained (except for the 
respective square-integrability requirements). The building blocks in (6.10) are the 
matrix-valued wavelets parameterized by E, whereas those in (6.11) are the matrix- 
valued plane-wave solutions e ip ' x II(p) parameterized by C. 

7. Interpretation of the Wavelet Parameters 

Our goal in this section is twofold: (a) Reduce the wavelets * 2 to a sufficiently simple 
form that they can actually be visualized, and (b) use the ensuing picture to give a 
complete physical and geometric interpretation of the eight complex space-time param- 
eters z G T labeling ^ z . That the wavelets can be visualized at all is quite remarkable, 
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since *& z (x') is a complex matrix- valued function of x' G R 4 and z G T. However, 
the symmetries of Maxwell's equations can be used to reduce the number of effective 
variables one by one, until all that remains is a single complex-valued function of two 
real variables, whose real and imaginary parts can be graphed separately. 

We begin by showing that the parameters zeT can be eliminated entirely. Recall 
that * z (x') is the boundary value of the reproducing kernel, according to (5.11) and 
(3.7): 

ty z (x') = J lim [K(x' + iey' | z) + K(x' - iey' | z)} 

= I lim [%' • y)L(x' + iey' - z) + 6(-y' ■ y)L(x' - iey' - z)] (7.1) 

= \ W ■ v) + 0(-y' ■ y)] Mx' -z) = \ Mx' - z). 

Hence 

ty x+iy (x') = l - L(x' -x + iy) = y iy (x' - x), (7.2) 

and *& x +iy is a translated version of . It therefore suffices to examine only with 
y G V . Eq. (4.7), combined with Tl(—p) = Tl(p), shows that *& z (x')* = *&z(x'), hence 
it suffices to look only at y G V\_. To reduce the number of parameters still further, 
we use the fact that Maxwell's equations are invariant under Lorentz transformations, 
and this invariance implies certain transformation properties for the wavelets. The 
covariance of the wavelets under the Lorentz group and, more generally, under the 
conformal group, will be studied in detail elsewhere. Here we remark only that Lorentz 
transformations relate all the wavelets with equal values of y 2 , hence it suffices to study 
only with y = (0, s) and s > 0. The physical significance of this will be discussed 
below. Finally, note that T(ap) = T(p) for any a > 0, since v(ap) = ap/ap = v(p). 
Hence U(ap) = II(p), and (4.7) implies that 

*o, ls (x')=s~ 4 *oAx'/s)- (7-3) 

Thus all the wavelets z G T, can be obtained by space-time translations, Lorentz 
transformations and scalings from the single "mother wavelet" 

&{x) = & ,i{x) = 2 f dp pi e~ P0 e ip ' x U(p). (7.4) 
Jc+ 

(Of course, any one of the \I/ 2 's can equally be chosen as the mother!) In particular, 
the wavelets parameterized by (x, is) G E are 
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* x , ls (x',t')=«- 4 *(^ 7 ^,Q. (7.5) 

Let [*(x, t)] mn denote the matrix elements of *(x, t). By (5.10), with wo = t+i, w = x 
and r = |x|, we have 



r _/ ,xi 1 <5 mn [(t + z) 2 + r 2 ] - 2x m x n + 2i(t + i) Y^k=i 

t)\ mn - — _____ . u . bJ 

This is still a complex matrix- valued function in R 4 , hence impossible to visualize 
directly. We now eliminate the polarization degrees of freedom. Returning to the 
Fourier representation of solutions, note that if f(p) already satisfies the constraint 
(2.12), then U(p) f (p) = f (p) and (4.10) reduces to 

[ d 3 xcfe|F(x,zs)| 2 = [ ^-|f( p )| 2 = ||F|| 2 . (7.7) 

J E JC Po 

Define the scalar wavelets by 

* z (x')= [ dp 2p 2 6(p-y)e- p < x '-^ (7.8) 
Jc 

and the corresponding scalar kernel K : T x T — > C by 

K{z' \z)= [ dp 4p 2 e(p ■ y') 6{p ■ y) e ip < z '~ s \ (7.9) 
Jc 

Then (7.7), with essentially the same argument as in the proof of Theorem 1, now gives 
the relations 



F(z') = 


/ d 3 x 
Je 


ds K(z' x, — is) F(x, is) pointwise in z £ T, 


F(x') = 


[ rf 3 x 
Je 


ds ^ Xjis (x')F(x,is) a.e. 


F = 


[ rf 3 x 
Je 


ds ^x^s F(x, is) weakly in H. 



(7.10) 



The first equation states that K(z' \ z) is still a reproducing kernel on the range W of 
Re : H — > C 2 {E). The second and third equations state that an arbitrary solution 
F G 7i can be represented as a superposition of the scalar wavelets, with F = Re F as a 
(vector) coefficient function. Thus, when dealing with coefficient functions in the range 
W of Re, it is unnecessary to use the matrix- valued wavelets. The main advantage of 
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the latter (and a very important one) is that they can be used even when the ceofhcient 
function 3>(x, is) does not belong to the range W of Re, since they project to W. 

The scalar wavelets and kernel were introduced and studied in Kaiser [12, 13]. They 
cannot, of course, be solutions of Maxwell's equations precisely because they are scalars. 
But they do satisfy the wave equation, since every component of \I/ 2 does so. To see 
their relation to the corresponding matrix quantities, note that Tl(p) is a projection 
operator of rank 1, hence Trace U(p) = 1. Taking the trace on both sides of Eqs. (4.7) 
and (4.20) therefore gives 

V z (x') = Trace V z (x'), K{z' \ z) = Trace K(z' | z). (7.11) 

Taking the trace amounts, roughly, to averaging over polarizations. The trace of the 
mother wavelet \I/ is 

Trace *(x) = -L = *M), r=|x|. (7.12) 

Because it is spherically symmetric, \l/(r, t) can be easily plotted. Its real and imaginary 
parts are shown in Figs. 1 and 2. These figures confirm that \1/ is a spherical wave 
converging towards the origin as — oo < t < 0, becoming localized in a sphere of radius 
v^3 around the origin at t = 0, and then diverging away from the origin as < t < oo. 
Figure 3 shows \I/(r, 0), which is real. Even though \l/(r, 0) does not have compact 
support (it decays as r -4 ), it is seen to be very well localized in |x| < \/3. 

Now that we have a reasonabe interpretation of ^f x +iy with y = (0, s), let us 
return to interpret the wavelets with y^0. Suppose y = (y, s) £ V\_ , and let v = y/s. 
Then y e V+ implies |v| < 1. (We have chosen units of length and time in which 
the speed of light c = 1; for general units, |v| < c.) Hence we can perform a Lorentz 
transformation to a reference frame moving with velocity — v relative to the original 
frame. In the new frame, y has coordinates y' = (0, -\J s 2 — |y| 2 ), hence our wavelet 
has a stationary center. Returning to the original frame, we conclude that ^ x +i y is a 
wave whose center is moving with the uniform velocity v = y/s. It is a Doppler-shifted 
version of the stationary wavelet with y' = 0. Thus each of the eight real parameters 
z = (x, t) + i(y,s) G T has a physical and geometric significance: x and t give the 
location and time at which \I/ 2 is localized; v = y/s gives the velocity of its center; 
\s\ gives its scale (width at time t), and the sign of s gives its helicity Since all these 
parameters, as well as the wavelets which they label, were a direct consequence of the 
extension of the electromagnetic field to complex space-time, it would appear that T is 
a natural arena in which to study electrodynamics. 
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8. Moving and Accelerating Wavelet Representations 

The construction of the electromagnetic wavelets has been completely unique, in the 
following sense: (a) The inner product (2.31) on solutions is uniquely determined, up 
to a constant factor, by the requirement that it be Lorentz-invariant. (b) The analytic 
extensions of the positive- and negative-frequency parts of F to 7+ and 71 , respectively, 
are certainly unique, hence so is F(z). (c) The evaluation maps £ Z F = F(z) are unique, 
hence so are their adjoint * z = £*. On the other hand, the choice of the Euclidean 
space-time region E as a parameter space for expanding solutions is rather arbitrary. E 
may be regarded as the group of space translations and scalings, acting on real space- 
time by g Xt i s x' = sx' + (x, 0). As such, it is a subgroup of the conformal group C, which 
consists of space-time translations, scalings, space rotations, Lorentz transformations 
and special conformal transformations. E is invariant under space rotations but not 
under time translations, Lorentz transformations or special conformal transformations. 
This non-invariance can be exploited by applying any of the latter transformations to 
the resolution of unity (4.16) and using the transformation properties of the wavlelets. 
The general idea is that when g e C is applied to (4.16), then another such resolution of 
unity is obtained in which E is replaced by its image gE under g. If gE = E, nothing 
new results. If g is a time translation, then the wavelets parameterized by gE are all 
localized at some time t ^ rather than t = 0. If g is a Lorentz transformation, then 
all wavelets with z G gE have centers which move with a uniform non-zero velocity 
rather than being stationary. Finally, if g is a special conformal transformation, then 
gE is a curved submanifold of T and the wavelets parameterized by gE have centers 
with varying velocities. This is consistent with results obtained by Page [20] and Hill 
[9], who showed that special conformal transformations can be interpreted as mapping 
to an accelerating reference frame. 

As a possible application, consider an electromagnetic pulse reflected or emitted 
by a moving object. After the reflection time, and far away from boundaries, the pulse 
may be approximated by a solution of Maxwell's equations in free space, hence it can be 
analyzed as in Section 4 using wavelets with stationary centers. However, the analysis 
is likely to be more efficient (i.e., have fewer significant coefficients) if it is made in 
the reference frame in which the reflecting object is at rest, with the reflection time as 
the initial time of localization. From the viewpoint of the receiver, this means that a 
representation with "co-moving" wavelets should be used instead of one with stationary 
centers. The details will be presented elsewhere. 

I thank R. F. Streater for his hospitality at King's College, where we had some 
helpful discussions concerning the helicity of the electromagnetic wavelets. 
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